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Abstract 

We present a pion photoproduction model on the free nucleon based on an Effec- 
tive Lagrangian Approach (ELA) which includes the nucleon resonances (A(1232), 
N(1440), N(1520), N(1535), A(1620), N(1650), and A(1700)), in addition to Born 
and vector meson exchange terms. The model incorporates a new theoretical treat- 
ment of spin-3/2 resonances, first introduced by Pascalutsa, avoiding pathologies 
present in previous models. Other main features of the model are chiral symme- 
try, gauge invariance, and crossing symmetry. We use the model combined with 
modern optimization techniques to assess the parameters of the nucleon resonances 
on the basis of world data on electromagnetic multipoles. We present results for 
electromagnetic multipoles, differential cross sections, asymmetries, and total cross 
sections for all one pion photoproduction processes on free nucleons. We find overall 
agreement with data from threshold up to 1 GeV in laboratory frame. 
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1 Introduction 



In spite of the fact that Quantum Chromodynamics (QCD) is regarded as 
the theory of the strong interaction, in the energy regime of the mass of the 
nucleon and its resonances a perturbative approach is not suitable. Thus, we 
have to rely on an effective approach if we are interested on the properties 
of nucleon resonances and processes where they are involved - mainly meson 
production which is the dominant decay channel. This paper is devoted to 
pion photoproduction from the nucleon, a classical topic within nuclear and 
particle physics, which has been proved as one of the best mechanisms to 
study the nucleon and its resonances as well as to study the role of the pion 
and resonances in nuclei [1]. 

The experimental database [2,3,4] has been enormously increased thanks to the 
experiments carried out at LEGS (Brookhaven) [5,6] and MAMI (Mainz) [7] 
where photons are produced through laser backscattering and bremsstrahlung 
respectively. Because of this experimental effort, our knowledge of the A(1232) 
resonance region has been largely increased, though several discrepancies be- 
tween Mainz and Brookhaven analyses still remain [6]. Polarization observ- 
ables, differential cross sections and electromagnetic multipoles have been 
measured with a precision not possible a few years ago and a full descrip- 
tion of the amplitudes in the A(1232) region is now available. The database is 
expected to grow significantly once data from current experiments have been 
analyzed and when data from new laser backscattering facilities as GRAAL 
(Grenoble) and LEPS (Harima) become available. The last two facilities have 
started to run recently and operate at higher energies than LEGS. This situ- 
ation opens a lot of possibilities for research on nucleon resonances. 

In the last decades, pion photoproduction has been studied through many 
models and using various approaches to the description of the nucleon reso- 
nances. Among them there are, Breit-Wigner models [8,9], K-matrix [10,11], 
Effective Lagrangian Approach (ELA) [12,13,14,15], dynamical models [16,17,18,19], 
Breit-Wigner plus a Regge-pole type background to take into account the ex- 
change of heavier mesons [20], as well as quark models with pion treated as 
an elementary particle [21]. Although in one way or another all models are 
phenomenological, in this paper we adopt ELA method because we consider it 
appealing in many respects and it is the most suitable approach in the energy 
range from threshold up to 1 GeV in laboratory frame, where the main low- 
lying resonances live. This approach has proved to be a quite succesful tool to 
study pion photoproduction at low/threshold energy [22,23,24] and provides 
the most natural framework to extend the model to pion electroproduction 
[12], electromagnetic pion production in composite nuclei [25] and halo nu- 
clei [26], two pion photoproduction [27], meson exchange currents [28], and 
exclusive X(7,A'"7r)Y processes. 
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In the last years the Lagrangian description of spin-3/2 resonances has been 
greatly improved and many pathologies related to the pion-nucleon-resonance 
and 7-nucleon-resonance interactions have been overcome [29]. This fact, com- 
bined with the substantial enlargement of the pion photoproduction database, 
demands to revisit the topic and to make the most of these advances in or- 
der to improve our knowledge on nucleon resonances and 7-nucleon-resonance 
vertices as well as on the pion photoproduction process itself. 

We focus on the analysis of pion photoproduction process on free nucleons with 
the aim of establishing a reliable set of coupling constants and achieving an 
accurate knowledge on nucleon resonances. The latter are needed for further 
studies of resonances in nuclear medium as well as to study the structure of the 
nucleon through its excitations. This requires to develop a pion photoproduc- 
tion model and to study the parameters of the nucleon resonances within the 
model for further implementation in the calculations previously mentioned. 
In this regard, we consider this paper as a first step towards a deeper un- 
derstanding of the role of the pion and the resonances in more complicated 
processes. Our model is an improvement of the one in Ref. [12] where we have 
changed the spin-3/2 Lagrangians and explored other variations which allow 
us to achieve crossing symmetry and a better description of the resonance 
widths. The elements included in this model are nucleons, pions, photons, p 
and uo mesons, as well as all four stars status spin 1/2 and 3/2 nucleon res- 
onances up to 1.7 GeV in Particle Data Group (PDG, in what follows) [30]. 
Spin-5/2 resonances are not expected to play an important role in the data 
analysis carried out in this paper and are left to future exploration. 

The paper is organised as follows: in section 2 we provide the basic features 
such as conventions and normalizations for cross sections and amplitudes 
which will be used throughout the article. In section 3 we describe the full 
model and its features in detail, stressing crossing symmetry and the spin-3/2 
treatment which avoids well-known pathologies of previous models. In sec- 
tion 4 we show results for multipoles, differential cross sections and remaining 
physical observables. We also provide the reader with all the parameters of the 
model explaining how they have been determined. In section 5 we summarize 
the main conclusions and results. 



2 Kinematics, cross section and amplitude decomposition 



Notation for kinematics is set to A; = (i?^, k) for the outgoing pion, q = {E^, q) 
for the incoming photon, p = {E,p) for the incoming nucleon, and p' = {E',p') 
for the outgoing nucleon. Mandelstam variables are defined as usual [31] 
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u^(p' -qf = (p-kf , 
t^{k-qf^ip-p'f . 



(1) 
(2) 
(3) 



The photon polarization vectors in spherical basis are 

^A,=±i = T;^(0,l,±^,0) . (4) 

Following conventions and normalization of Rcf. [31], the differential cross 
section can be written in the center of mass (cm.) reference system as 

Whenever a kinematical quantity appears starred it is defined in the cm. ref- 
erence frame. In particular the cm. absolute values of the photon and the pion 

J -* I I ~* I 

momenta are denoted by q* and A;*, which stand for |g*| and \k*\ respectively. 
The transition probability is 

W?=\Y. (6) 

where .4aiA2A^ is the invariant amplitude, with photon polarization A^, initial 
nucleon helicity Ai, and final nucleon helicity A2. We use the following isospin 
decomposition 

A = x\ (^S + ^"^[r,-, ra] + vr.Xi , (7) 

where for simplicity we have dropped helicity subindices (another two isospin 
decompositions are used and are introduced in appendices). 

The isospin decomposition can be related to the physical amplitudes 



A {jp ^ p7r°) =A+ + A\ 
A{yn^n7r'')^A+ - A\ 
A{yn^pn-)^V2(^A'^ - A-) , 
A (tp ^ n7r+) ^V2[A^ + A') . 



(9) 
(10) 

(11) 
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(a) (6) (c) id) 



Fig. 1. Feynman diagrams for Born terms: (a) direct or s-channel, (6) crossed or 
■u-channel, (c) pion in flight or t-channel, and {d) Kroll-Rudermann (contact). 

For completeness we specify the conventions adopted throughout this arti- 
cle. Metric tensor: g'^'^ = diag{l, —1, —1, —1); Levi-Civita tensor: 60123 = 1, 
^123 = 1; Pauli matrices: [Tj,rfc] = 2ejkiTi] Dirac-Pauli matrices: {7^,7*^} = 
2g'"', Y"" = \ [7^,7l, 7'"'^" = \ {-i^'YT - TYl^). 75 = ^7o7i7273; Electro- 
magnetic field: F'^'^ = d^Ay - d'' M , F'^'^ = ie^^^^^Fa^, where is the photon 
field. 



3 The model 



In this section we present a complete description of the model and its features. 
Using as starting point Weinberg's theorem [32], we construct a fully relativis- 
tic, chiral symmetric, gauge invariant, and crossing symmetric model based on 
suitable effective Lagrangians for particle couplings. From these Lagrangians 
we obtain the invariant amplitudes and physical observables. This procedure 
has been adopted in many papers, i.e. [12,13,14,15], and has been proved to be 
a succesful way to treat the pion photoproduction process. However, previous 
works had pathologies in the description of the spin-3/2 particles which are not 
present in our model. The basic idea is to build consistently the most general 
Lagrangians for vertices, taking into account all possible symmetries (crossing 
symmetry, gauge invariance, chiral symmetry), and to use Feynman rules to 
obtain invariant amplitudes which can be related to physical observables. The 
model can be split into three different types of contributions: Born terms (Fig. 
1), vector mesons exchange (Fig. 2, diagram (e)), and spin-1/2 and spin-3/2 
nucleon resonance excitations (Fig. 2, diagrams (/) and i^g)). There is no con- 
tribution from a meson exchange because of charge conjugation violation of 
the (T7r7 coupling [11]. We consider that all the relevant degrees of freedom 
are taken into account except perhaps spin-5/2 resonances. Our choice of La- 
grangians is explained and justified in the forthcoming subsections. All the 
invariant amplitudes can be found in appendix A. 



3.1 Born terms 



Born terms are the Feynman diagrams shown in Fig. 1 in which only pi- 
ons, photons, and nucleons are involved. We start from the free Lagrangians 
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for pions (Klein- Gordon) and miclcons (Dirac) and a phenomenological pion- 
nucleon interaction. This last interaction is chosen as a pseudovector (PV) 
coupling to the pion because it is the lowest order in derivatives compatible 
with the low energy behavior of the pion and chiral symmetry [22,33] 

C^NN = — A^7M75r,- (a'^TT,) TV , (12) 

where is the mass of the pion, /^jv is the pseudovector coupling con- 
stant, and the sign is fixed by phenomenology. According to [34] /ttAt is set 
to fl^jAn = 0.0749. The use of the PV coupling for the pion in our effective 
Lagrangian grants that low energy theorems of current algebra and partially 
conserved axial-vector current (PCAC) hypothesis are incorporated in the 
model. 

The electromagnetic field is included in the usual way by minimal coupling 
to the photon field (9" —>■ + ieQA^] where Q is the charge operator) and 
taking into account phenomenologically the anomalous magnetic moment of 
the nucleon 



^ = -m'^^^l {^'"^ + 7«/.iVF«^ . (13) 

¥2^^ is defined as the ratio between isospin-scalar and isospin-vector form 
factors {F^'^ = Fi/F^). 

The interacting Lagrangian for Born terms is 

Cbotu = - ieF,i"e,-fe37r,- {d^T^k) - eA'^F^N^J^ (f^"" + T3) N 

- ieF^^A^^N^^^.^Tj, n]7rjN (14) 

where e is the absolute value of the electron charge, is the pion form factor 
and Fy — Fj — F", F^ = Fj -|- Fj* are the isovector and isoscalar nucleon 
form factors which at the photon point (g^ = 0) take the values F^ = FY = 1, 
= kP + k"" = -0.12, F^ = kP - k"" = 3.70. We set F^ = F^ in order to 
ensure gauge invariance. It is straightforward to check gauge invariance of the 
amplitudes in appendix A performing the replacement A^^ — > q^^. 
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(e) (/) {g) 

Fig. 2. Feynman diagrams for vector meson exchange (e) and resonance excitations: 
(/) direct or s-channel and {g) crossed or n-channel. 

3.2 Vector mesons 



The main contribution of mesons to pion photoproduction is given by p 
(isospin-1 spin-1) and cu (isospin-0 spin-1) exchange. The phenomenological 
Lagrangians which describe vector mesons are [9,12] 



^uj — —Fu,nnN 



Lp — —FpjsjjsjN 



'2M 
2M 



7« - ^TTir^ic^P^^ 



^-N + ^^F^, (d'^nj) Sj^cu^ , (15) 
r,p«Ar + ^F,.(5%.)p;. (16) 



Often the tt^V couphng is written as £ = ^^^F^i^V'^'Tr where = d/^V - 
d^V^ and = ff'^u'^ [9]. Both couphngs yield the same amphtude. 



3.3 Spin- 1/2 nucleon resonances 



In the model we deal with three different kinds of resonances with spin-|: Sn, 
S31, and Pii; and we need Lagrangians and amplitudes to describe their be- 
havior. The most simple isobar is isospin-| spin-| (Sn) which can be described 
by the following Lagrangian 

h — iG — 

>Csn = -^iV7,r,Ar*a-7r,- - jjjN'y^f.'y, {gs + gv^) N*F^^ + HC , (17) 



where HC stands for hermitian conjugate, h is the strong coupling constant 
which can be related to the width of the resonance decay into a nucleon and 
a pion, and f^^ = 92.3 MeV is the leptonic decay constant of the pion. gy and 
gs stand for the resonance isovector and isoscalar form factors respectively. 
They are defined as gv — Qp — Qn and gs — Qp + Qn, where subscripts p and 
n stand for the resonances originating from the proton and the neutron, and 
can be related to experimental helicity amplitudes at the photon point as will 
be seen in the next sections. The pion coupling has been chosen pseudovector 
in order to obtain the right low energy behavior and consistency with Born 
terms. The coupling to the photon used preserves gauge invariance. 
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The next isobar is isospin-| spin-| (S31). Because of isospin we need to define 
isospinors as in Ref. [35] 



N* 



N* 



/ 



^N*+ - N*- 
I 




(18) 
(19) 
(20) 



In this basis, and under the same conditions as those for previous isobar, the 
S31 Lagrangian is 



N^ap^,N;F'''' + HC . 



(21) 



Just one electromagnetic coupUng constant is needed here because only the 
isovector part of the photon couples to the nucleon to produce an isospin- 1 
field. 

The Pii Lagrangian is closely related to Sn being parity the main change. 
This change is due to the angular momentum of the resonance, which implies 
a change in the parity of the coupling 

h _ op _ 

Cp,, = --7V7«75T,JV*a"7r,- + ^%a/3 (^5 + gvrs) N*F''^ + HC . (22) 



3.4 Spin-3/2 nucleon resonances 



The choice of spin-3/2 nucleon-resonance couplings is one of the main im- 
provements of the present model compared to former ones. The choice that 
we use here is motivated by previous studies that identified pathologies in 
former spin-3/2 couplings. In what follows we provide a detailed comparison 
of both traditional (off-shell extension) and gauge invariant (GI) couplings, 
which exhibits the virtues of the choice adopted here. With regards to the 
traditional coupling, we restrict the discussion to the P33 (A) resonance and 
its coupling to the pion and the nucleon, although it affects similarly to other 
spin-3/2 resonances. 
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3.4-1 Traditional A-nucleon-pion coupling 



The basis of the traditional point of view is the seminal paper by Nath, 
Etemadi, and Kimel [36], based on the articles by Peccei [35] in the late six- 
ties which dealt with this coupling. Peccei worked out a chiral Lagrangian 
with a pseudovector coupling to the pion, to ensure the low energy behavior, 
based upon the invariancc of the A free field under the point transformation 
— > A^ — jji^jpA'^ and the ansatz ■y'^O^u = 0. Given the most general 
Lagrangian C — hA'^O^^Nd'^-Kj we obtain the well known Peccei Lagrangian 
[35] 

Cpeccei = ihA'^ {^9aP - lalp) Nd^TTj + HC . (23) 

Restrictions such as Peccei's ansatz are needed in order to reduce the number 
of degrees of freedom (DOF) of the spin-3/2 field. When a massive spin-3/2 
particle is described within the Bargmann and Wigner equations [37] , a prob- 
lem of extra DOF arises because a vector-spinor has sixteen componentes 
whilst only four are needed. These constraints naturally emerge in the free 
theory thanks to the Euler-Lagrange [36] or the Hamiltonian formalism [29], 
but for interacting particles the picture is not so straightforward and addi- 
tional restrictions have to be imposed. 

Nath et al. [36] proved that Peccei's ansatz was too restrictive, developing 
a generalization which - despite of its many pathologies [29,36,38,39] - has 
become the traditional and most popular approach to interacting spin-3/2 
particles for the last thirty years. 

The starting point of Nath et al. is the massive spin-3/2 free theory, which 
can be found in Refs. [35,36,38]. The following Lagrangian is defined 

(i9^7'' - M*) g^p + iu {-fadp + -y^da) 

n (24) 

+ ^ (3a;2 + 2uj + l) 7a5^7M7/3 + M* {Scu^ + Scu + l) 7^7^] A'^ , 
where a; 7^ — | and the Lagrangian is invariant under the point transformation 



A'^ ^ A'^ + a7''7"A, , (25) 
uj — 2a , , 

with a 7^ — |. Subsidiary constraints ^^Ai^ — and 9^ A'* = appear in order 
to reduce the number of DOF to four, as expected for a spin-3/2 particle. 

A detailed description of the DOF counting technique is given in reference 
[40]. The parameter u does not affect physical quantities, so that one is free 
to set it to the most convenient value, usually u = —1 which recovers the 
Rarita-Schwinger theory [41]. 
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The point transformation of Eq. (25) does not affect the spin-3/2 content 
of the free field because of the constraint 7^A'' = 0, but for interacting A 
particles this constraint does not apply, and the excess of DOF shows up as a 
contribution to the spin- 1/2 sector. 

The most general interacting Lagrangian containing only first-order derivatives 
of the pion field and consistent with (24), (25), and (26) is given by 

Lint = xA« (^„/3 + a^oHp) Nd^n + HC , (27) 

where x is a coupling constant and a is called the off-shell parameter, which 
can be set to different values. This is named off-shell extension framework. If 
a = — I we recover Peccei theory. This family of Lagrangians has been widely 
used in pion-nucleon scattering [15,42], pion photoproduction [12,13,14,15,16,17,35] 
and compton scattering [15,43] in the A-region, as well as for the description 
of meson exchange currents [28]. The off-shell parameter can be set to a fixed 
value, a = — 1 [36], a = —\ [12,35] or just let it run freely [13,43] in order to 
get the best possible fit. 

However, it is not possible to remove the spin- 1/2 sector from the amplitude for 
any value of a [38] . The physical meaning of the off-shell parameter is unclear 
and could be considered just as a free parameter with a fuzzy physical meaning 
set only for fitting improvement. The disadvantage is that there is a heavy 
dependence of the coupling constants on the off-shell parameter, as was proved 
by Feuster and Mosel [13]. Other pathologies related to Eq. (27) coupling are: 
quantization anomalies (except for a = —1), so that the naive Feynman rules 
we read from the Lagrangian are no longer valid [29,36]; Johnson-Sudarshan 
(JS) problem (nonpositive definite commutators) [44,45] and Velo-Zwanziger 
(VZ) problem (acausal propagations) [45,46]. 

A consistent theory for interacting spin-3/2 particles is expected to be free 
of such problems. This theory has been developed in recent years and will be 
detailed in the next paragraphs. 

3.4-2 Gauge invariant couplings 

A different approach to massless fields of arbitrary spin A was developed in the 
seventies. It was proved that the massless theory obtained from the massive 
one has a simple structure for both integer [47] and half-integer [48] spin fields, 
even if the massive theory is rather complicated. The free massless Lagrangians 
for half-integer spin fields can be obtained just from first principles requiring 
the action to be invariant under the gauge transformation ip ip + Sip, where 
Sip = drj [49,50], ip is a tensor-spinor with rank i which stands for the parti- 
cle and 7] a complex tensor-spinor field with rank i — 1. For a spin-3/2 field 
Sip^ = dfj,r], with ipn a vector spinor and rj a spinor field. This gauge condition 
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reduces the number of DOF of the spin-A field to 2 - hehcity states —A and 
+A - as it is required for a massless particle. In this framework, it is quite sim- 
ple to build consistent interactions for half-integer spin fields as suggested by 
Weinberg and Witten [51] just enforcing them to fulfill this gauge invariance 
condition. For example, the spin-3/2 ■^^ field should appear in the interaction 
as "the spin-5/2 as dij,d,,iljpa-dij,daippa-dpd,yiljf,a+dpdaipf,u, and, 

more generally, an arbitrary spin-A tensor-spinor field as the antisymctriza- 
tion of daida2 ■ ■ ■ ^a(;^_i/2)V'/3i/32 - /3(A-i/2)- Thus the vertices 0^^^^ of the Feynman 
diagrams for massless spin-3/2 particles will fulfill the condition p^C'* " = 
where p is the four- momentum of the spin-3/2 particle, /i the vertex index 
which couples to the spin-3/2 field, and the dots stand for other possible in- 
dices. This is what is called GI coupling scheme. 

We apply this procedure to the A case. We start from Lagrangian (24) for a 
free massless spin-3/2 particle. For a; = — 1 it can be written as 

^3/2,massless = i^nY^^dai^u ■ (28) 

The inclusion of the mass term 

^3/2,massive = ^3/2,massless — M^ll) ^^y^^ ifjy , (29) 

breaks gauge symmetry, raising the number of DOF from 2 to 4 as it should 
be. 



Let us now consider the interaction. For an interacting massless spin-3/2 par- 
ticle we write the Lagrangian 

'~'3/2,massles + Ant ■ (30) 

The interaction has been built within the GI coupling scheme and can be 
written as [29] 

Ant = ^iJ" + HC , (31) 
where has no dependence on and gauge invariance imposes d^J^ — 0. 
The inclusion of the mass term - if it is properly done as in (29) - breaks 
gauge symmetry increasing the number of DOF of the spin-3/2 field from 
2 to 4 and does not affect Cint [29,50]. Hence the number of DOF in the 
interacting massive field is the right one and no unphysical components are 
present. Focusing on our photoproduction model, we are interested in two 
couplings: the A to the pion and the nucleon, and the A to the photon and 
the nucleon. The simplest consistent AA'"7r-couphng is [29] 

Ant = -j^Ne,^xpl^^' (^'^iVf ) {d%) + HC . (32) 

We have to clarify that the vector coupling to the pion is a consequence of GI 
prescription. Whithin this prescription, the scalar coupling to the pion gives 
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no contribution to the amplitude [29]. 



Concerning tiie AN7 coupling, Jones and Scadron [52] suggestion has been 
widely used in the {Gx,G2) decomposition with [11,12,13,14,43] or without 
[16,17] off-shell extension. Another decomposition {Ge, Gm), based upon the 
same idea as the Sachs form factors for the nucleon [53], is also possible. 
This decomposition is directly connected to physical quantities, as electric 
and magnetic multipoles, in particular to the E2/M1 ratio which is of great 
interest from both experimental and theoretical points of view [6,19,24]. This 
second decomposition is consistent with the GI approach and can be written 
as [54] 



C = N 

2M{M + M*) 



WiF^u + g2l'F^,\ (a'^TV-) + HC , (33) 
where gi and g2 can be easily related to Ge and Gm [40] by 



Ge — —- 



IM* -M 
2M* + M^^ ' 



_ 1M*-M 



(34) 

(35) 



Other possible consistent choices can be found in Refs. [40,55]. 

GI couplings have been proved to be free of the pathologies which are inherent 
to the traditional scheme. No anomalies are found in the quantization; neither 
JS nor VZ problems appear; and no spin-1/2 sector arises when the invariant 
amplitudes arc calculated [29]. Moreover, Pascalutsa and Timmermans [40] 
claim that DOF counting is the reason why GI couplings are consistent while 
the off-shell extension couplings of Nath et al. are not. They blame the un- 
physical extra components for the appearance of pathologies. Both, GI (32) 
and traditional (27) couplings, provide the same result on-shell (if wc set prop- 
erly the coupling constants). However, their off-shell behavior is completely 
different. 



Based on the previous discussion, the P33 Lagrangian that will be used in this 
work is 



h 



'Pas 



+ 



3e 



2M (M + M*) 
+ HG . 



N 



(36) 



Prom this Lagrangian it is straightforward to obtain phenomenological La- 
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grangians for other spin-3/2 resonances. To obtain the P13 resonance La- 
grangian from (36) only an isospin change is needed 



J 



r,N*-, J = 1,2,3 



(37) 



for the strong vertex, and 



N* 



9j 



1,2; 



(38) 



for the photon vertex. 
Thus the Lagrangian is 



h 



+ 



(9t + 9Xr:i)F^u 



4:M{M + M*) 
+ {9I + 92 rs) I'F, 



{d^N*") + HC 



(39) 



Lagrangians for D33 and D13 resonances are obtained easily from P33 and P13. 
We only need to change the parity of the coupling placing an overall 75 



-^033 — 



h 



+ 



3e 



2M (M + M*) 
+ HC, 



N 



(40) 



>^Dl3 — 



h 



+ 



3e 



AM (M + M* 
+ {92 + 92 rs) F^, 



N 



i [9i + 9irs) F^^-f5 



(41) 



{d^N*") + HC . 



Although we restrict ourselves to spin-3/2, it is clear that higher spin interac- 
tions can be built within the same theoretical framework. This is left to future 
works. 
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3.5 Propagators and widths 



With regards to the propagators of the resonances, for a spin- 1/2 resonance 
we use 



\^ - M*^ +iM*T{s,u)' 
and for the spin-3/2 propagator we use the Rarita-Schwinger propagator 



(42) 



X 



y2 - M*^ + iM*r{s,u) 
1 2 

-9aP + g7a7/3 + 



3M' 



(43) 



where v is the resonance four-momentum. A phenomenological width V (s, u) 
is included in the propagator denominator consistently with what is obtained 
if we dress it with pions [54,56]. 



The energy dependence of the width is chosen phenomenologically as 



(44) 



where j — tt, tttt, r] stands for the different decay channels and 

X, (s, u) ^ X, (s) + X, {u) - X, {s) X, {u) , 
with XjiJ) given by 



X, (/) 



\H 

l^jol 



2L+1 



2L+3 



e (/ - (M + m,)') , 



(45) 



(46) 



where L is the angular momentum of the resonance, © is the Heaviside step 
function, and 



\k^\ = - mf) - Am^M^/ (2V1) , (47) 

with mTTTT = 2m7r and \kjQ\ — \kj\ when I — M*^. 

This parametrization has been built in order to fulfill the following conditions 

(i) r = To at = M*, 

(ii) r — > when \kj\ — > 0, 

(iii) a correct angular momentum barrier at threshold \kj\'^^~^^, 

(iv) crossing symmetry. 
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This paramctrization of the width is an improvement over the one used in Ref. 
[12] and includes decays to rj and 27r which take into account inelastic channels 
[9] and condition (iv). The width contributes to both s and w-channels, so that 
crossing symmetry is preserved due to Eq. (45). In [13] the authors made an 
analysis of the energy dependence of the width. It was concluded that, as 
long as it provides a decrease of the width beyond the resonance position, the 
specific way in which Xj is parametrized is not so important. 



3. 6 Form factors 



For the numerical calculations we include form factors for Born terms and 
vector mesons, in order to regularize the high energy behavior of these terms. 
We choose form factors as suggested by Davidson and Workman [57] that 
allow to fulfill gauge invariance and crossing symmetry. Actually, Xj[s,u) in 
Eq. (45) also follows this choice. Thus for Born terms 

Fb{s, u, t) =Fi(s) + F^iu) + Fs{t) - Fi{s)F2{u) 

- F,(s)Fs(t) - F2(u)Fs(t) + F,(s)F2(u)Fs(t) , 

where. 





1+ 


[s-M-)'/K% 


F2{u) = 


1 + 






1 + 





1 -1 



(49) 
(50) 
(51) 



For vector mesons we adopt Fv{t) — F^{t) with the changes m^r — > my and 
As — > Ay. In order to have as few free parameters as possible in the numerical 
calculations we use the same A = Ag = Ky for both vector mesons and Born 
terms. For the resonance-pion-nucleon vertex, the form factor (s, u) has 
to be used for consistency with the width used in the propagator discussed 
previously. 

Models like the ones by Garcilazo and Moya de Guerra [12] and Feuster and 
Mosel [13] needed a cutoff in the w-channels of spin-3/2 resonances to obtain 
a good description of observables. This cutoff was needed because the high- 
energy contributions of these diagrams are not reduced by the denominator 
of the propagator. In Ref. [12] it was argued that the need of this cutoff could 
be justified by the two possible interpretations of the resonance excitation. 
From an effective field theory point of view, w-channels should be introduced 
with their full strength. On the other hand if we consider resonances as pure 
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t:N rescattering states (Chew-Low description), the u-channels contributions 
should be dropped. Hence, the cutoff was interpreted in Ref. [12] as a way to 
have an interplay between both descriptions. However, in this way crossing 
symmetry was broken. Our present model relies entirely on effective field the- 
ory, and we preserve crossing symmetry and there is no need for that cutoff 
in the w-channel amplitudes. Ought to the vector couplig to the resonance the 
M-channel amplitudes are suppresed by themselves which is a strong point in 
favor of the GI coupling. 



4 Results 



4.1 Study of the parameters of the model 

The first choice that has to be made is the nucleon resonances to be taken 
into account. We have included seven resonances: A(1232), N(1440), N(1520), 
N(1535), A(1620), N(1650), and A(1700) which are all the four star nucleon 
resonances in PDG up to 1.7 GeV and up to spin-3/2. Among four star reso- 
nances only spin-5/2 N(1675) and N(1680) resonances are left aside for future 
work. 

In a Lagrangian model, the determination of the parameters of a single res- 
onance is affected by the determination of the parameters of the other reso- 
nances. Thus, we have decided not to include three star resonances because 
their contribution would be very small and would introduce a sort of noise in 
the determination of the parameters. 

There are quite a number of parameters to be set in the model. Some of them 
are well known and established independently of the photoproduction data, 
such as nucleon and pion masses (M = 938.9175 MeV, rriT^o = 134.9766 MeV, 
m^ri = 139.5673 MeV), but some others have to be established from fits to 
the pion photoproduction data, namely electromagnetic coupling constants. 
In the forthcoming paragraphs we give the values of every parameter of the 
model as well as the procedures employed to establish them. 



4.. 1.1 Vector meson coupling constants 

Vector meson contributions are characterized by eleven parameters: m^^, F^^nn, 
K^j, Gu;7r7, ^po, mp±, FpNN, Kp, GpO^^^, Gp±T^^, and cutoff A. Masses are given 
by PDG and the n^V couphngs are related to the decay widths F^^y of PDG 
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Table 1 

Masses, widths, and branching ratios from Refs. [30,60] and from the speed plot 
calculation (see text). Masses and widths in MeV. We have taken Vt^t^/V = 1 — 
r^/r - r^/r. Subscripts P, V and SP stand for PDG [30], Vrana [60] and Speed 
Plot respectively. PDG masses and widths arc mean values. 





A(1232) 


N(1440) 


N(1520) 


N(1535) 


A(1620) 


N(1650) 


A(1700) 




1210 


1365 


1510 


1505 


1607 


1660 


1660 




1217 


1383 


1504 


1525 


1607 


1663 


1726 




1211 


1372 


1516 


1540 


1608 


1664 


1641 




100 


210 


115 


170 


115 


160 


200 




96 


316 


112 


102 


148 


240 


118 


r 

SP 


98 


290 


48 


107 


141 


159 


955 


la. 
Tp 


1.00 


0.65 


0.55 


0.45 


0.25 


0.72 


0.15 


Tp 




0.00 


0.00 


0.51 




0.06 




Tp 


0.00 


0.35 


0.45 


0.04 


0.75 


0.22 


0.85 




1.00 


0.72 


0.63 


0.35 


0.45 


0.74 


0.05 






0.00 


0.00 


0.51 




0.06 




Tv 


0.00 


0.28 


0.37 


0.14 


0.55 


0.20 


0.95 



[30] through the equation 



y— >7r7 



967r ml 



my, 



(52) 



We take from PDG the following values: = 782.57 MeV, mpo = 768.5 
MeV, mp± = 766.5 MeV, r^o^^ = 0.121 MeV {Gpo^^ = 0.1161), Tp±^^ = 0.068 
MeV {Gp±^^ = 0.0906), and T^^^ = 0.70476 MeV {G^^^ = 0.2804). Thus, 
only five constants remain unknown. One of them is the cutoff A which will 
be discussed later. The four remaining constants are taken from the analysis of 
nucleon electromagnetic form factors by Mergell, MeiBner, and Drechsel [58]: 
Fpj^j^ = 2.6, Kp = 6.1 ± 0.2, F^NN = 20.86 ± 0.25, and = -0.16 ± 0.01, 
which compare well to the data, including the latest experiments at Jefferson 
Lab [59]. 



4-1.2 Masses and widths of the nucleon resonances 

We have used three different sets of masses and widths of the nucleon reso- 
nances (Table 1): First, the PDG values [30]; second, the multichannel analysis 
of Vrana et al. [60] ; and third, the speed plot (SP) calculation that we explain 
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A(1232) 



N(1535) 



N(1 650) 



1600 1620 1640 1660 1680 



Fig. 3. Speed plot examples. Left figure shows the A(1232) speed plot. Right figure 
shows the speed plot for the Sn region. Data have been taken from SAID database 
for -kN scattering [4]. 

below. For the partial decay widths we have two different sets, one from PDG 
and one from Vrana et al. which hes within the PDG error bars. The Vrana 
et aL set of partial decay widths has been chosen for the SP calculation. 



Masses and widths of nucleon resonances can be obtained from nN partial 
wave analysis using the speed plot technique [61]. First we define the speed 

by 

SP {W) = \dT {W) /dW\ , (53) 

with W — y/s and 

T{W)^^.[ei^p{2iS{W))-l] , (54) 

where T{W) is the dimensionless resonance partial wave amplitude and S {W) 
its phase. 

To calculate masses and widths we have taken phases from the current solution 
of the SAID ttN partial wave analysis [4]. In figure 3 we show SP{W) for P33 - 
in the A(1232) region - and Sn multipoles. The position of the peak provides 
the pole mass, and the height provides the width: H = 2/T. 

The baryon resonances show up clearly and the calculation is straightforward. 
The only problem is related to the existence of a background which induces 
a phase shift in the ttN phases. In the region of the peak this phase shift can 
be considered approximately constant and its effect in SP{W) is negligible. 

Our fits to photoproduction data shown in the next subsections are clasified 
according to six sets of parameters which are given in table 2. Sets 7^1 and 
are based on PDG values for masses and widths; set #2 and set #4 are based 
on Vrana et al. [60] ; and sets #3 and #6 are based on the SP calcualtion. 
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Table 2 

Specifications of the parameter sets. Masses, widtlis, and A are in GeV. Tlie coupling 
constants for the vector mesons are dimensionless. We provide also the x^/Xpdg 
order to compare fits. 



Set 


#1 


#2 


#3 


#4 


#5 


#6 


Masses & Widths 


PDG 


Vrana 


SP 


PDG 


Vrana 


SP 


^FSI 


Yes 


Yes 


Yes 


No 


No 


No 


X^/Xpnc 


1 


0.53 


0.60 


9.30 


5.57 


4.56 


A 


1.121 


1.050 


1.040 


1.494 


0.951 


0.962 




6.30 


6.30 


6.30 


6.30 


5.90 


5.90 




21.11 


21.11 


21.11 


20.61 


21.11 


21.11 




-0.17 


-0.17 


-0.17 


-0.15 


-0.17 


-0.17 


M* [A(1232)] 


1.209 


1.215 


1.209 


1.210 


1.215 


1.209 


r[A(1232)] 


0.102 


0.098 


0.100 


0.102 


0.094 


0.099 


M* [N(1440)] 


1.385 


1.381 


1.370 


1.385 


1.381 


1.370 


r [N(1440)] 


0.160 


0.318 


0.292 


0.260 


0.314 


0.288 


M* [N(1520)] 


1.505 


1.502 


1.514 


1.505 


1.502 


1.514 


r [N(1520)] 


0.110 


0.110 


0.050 


0.110 


0.110 


0.050 


M* [N(1535)] 


1.495 


1.527 


1.542 


1.495 


1.523 


1.538 


r [N(1535)] 


0.250 


0.104 


0.109 


0.099 


0.100 


0.109 


M* [A(1620)] 


1.590 


1.605 


1.606 


1.620 


1.605 


1.606 


r [A(1620)] 


0.100 


0.150 


0.143 


0.100 


0.150 


0.143 


M* [N(1650)] 


1.680 


1.665 


1.666 


1.640 


1.665 


1.666 


r [N(1650)] 


0.150 


0.238 


0.157 


0.150 


0.242 


0.157 


M* [A(1700)] 


1.620 


1.728 


1.639 


1.620 


1.728 


1.639 


r [A(1700)] 


0.250 


0.120 


0.957 


0.250 


0.120 


0.957 



The strong coupling constants (/I's) of the resonances are obtained from equa- 
tions in appendix B, using the partial decay widths into one pion of the res- 
onances. We choose all the strong coupling constants to be positive, thus the 
overall sign of the amplitude of each resonance depends on the sign of the 
electromagnetic coupling constants. 
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4.1.3 Electromagnetic coupling constants of the nucleon resonances 

At this point, only the electromagnetic coupling constants and the cutoff A 
remain undetermined. The best way to establish them is by fitting to pion 
photoproduction experimental data. Among all the observables (cross section, 
asymmetries, etc.) for pion photoproduction, the set of data we have chosen 
is the one given by the current SAID multipole energy independent solution 
[2,3,4]. There are two main reasons for this choice. First, electromagnetic mul- 
tipoles are directly related to the amplitudes and are more sensitive to coupling 
properties than are other observables. Deficiencies in the model show up much 
more clearly in multipolcs than in any other observable. Second, all the ob- 
servables can be expressed in terms of the multipoles, thus, if the multipoles 
are properly fitted by the model, so should be the other observables. The ex- 
plicit expressions for the multipoles in terms of the amplitudes can be found 
in appendix C. 

Another issue to take into account is unitarity. Models below the two pion 
production threshold fulfill Watson's theorem [62] to achieve unitarity using 
either t^N scattering phases [9], dynamical models [16,17,18,19] or K-matrix 
[11]. Beyond the two pion production limit, implementation of unitarity is un- 
clear and usually relies on experimental data and/or extensions of the methods 
applied below the two pion threshold. 

We would hke to note that, although our calculation seems to be at tree- 
level, it is not quite so due to the inclusion of the width and the form factors, 
which take into account higher order diagrams and structure effects. If we 
perform a truly tree-level calculation - straightforwardly from amplitudes of 
appendix A - we would find out that all the amplitudes are real and that it 
would be impossible to fulfill the unitarity condition SS'^ — 1, where S is the 
scattering matrix. In an effective Lagrangian perturbative model, unitarity 
should be restored by the inclusion of higher order diagrams. To avoid this 
tedious and difficult task, we adopt a phenomenological point of view. The 
main higher order effects can be taken into account including a width in the 
propagator, as we do in section 3.5 (which amounts to dress the propagator), 
and including also effective final state interactions (FSI). Once the width is 
included, unitarity restoration may be achieved through FSI. We can assume 
that it is possible to isolate the FSI effects factorizing the multipoles M. in 
the following way: 

where F^'^'^ is a phase factor that takes into account FSI, and i. stands for 
orbital angular momentum, 11 for parity, and / for isospin: 

(56) 



exp 



i5 



FSI 
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Then, the absolute value of the multipoles must be well reproduced by the 
model and only the phases of the multipoles remain unknown. We are inter- 
ested in the bare values of the coupling constants, so the best choice is to use 
directly the experimental phases. Hence, the multipole phase can be written 
as 

S''''"" = 5^i,th + S'/^Y , (57) 

where we call Swidth to the phase given by the calculated amplitudes and com- 
parison with experimental phase shifts {6^'^'^) provides us with the unknown 
final state interaction phase shifts SpgY- Phases S^'^'^^ are taken from the cur- 
rent energy dependent multipole solution of SAID analysis [2,3,4]. For each 
set of masses and widths we obtain two types of fits, one with and one without 
SAID phases. 

In order to fit the data and determine the best parameters of the resonances we 
have written a genetic algorithm combined with the E04FCF routine from NAG 
libraries [63] . Although genetic algorithms are computationally more expensive 
than other algorithms, in a minimization problem it is much less likely for them 
to get stuck at local minima than for other methods, namely gradient based 
minimization methods. Thus, in a multiparameter minimization like the one 
we face here it is probably the best possibility to search for the minimum. 
It is out of the scope of this paper to go through an explanation on genetic 
algorithms and details on them can be found elsewhere [64]. 

The function to minimize is the defined as 



where Ai^^^ stands for the current energy independent extraction of the mul- 
tipole analysis of SAID up to 1 GeV for Eq^, Mi_, Ei^, Mi_j_, E2-, and M2- 
multipoles in the three isospin channels / = ^ for the 7p — > ir^p pro- 
cess. AM^^P is the error and Al*'* is the multipole given by the model which 
depends on the parameters. These parameters are the cutoff A and the elec- 
tromagnetic coupling constants in Table 3, which are related to the helicity 
amplitudes A{ in Table 4 through equations given in appendix D. 

The minimization procedure applied is as follows: First the genetic algorithm 
has been run and when the convergence conditions were accomplished the 
E04FCF routine was used for fine tunning. The program has been run many 
times with different seeds in order to ensure that the minimum was not local. 
We have taken into account 763 data for the real part of the multipoles and 
the same amount for the imaginary part. Thus, 1526 data points have been 
used in the fits. 

In Tables 2, 3, and 4 we show results for the six different sets and provide the 
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Table 3 

Coupling constants of the resonances. The E2/M1 Ratio (EMR) of A(1232) is also 
given. All magnitudes are dimensionless. 









II 1 


#2 


#3 


#4 


#5 


#6 


A(1232) 


P33 


h 


0.764 


0.721 


0.757 


0.759 


0.706 


0.753 






91 


6.061 


5.574 


5.630 


6.254 


5.382 


4.984 






92 


2.414 


1.187 


1.123 


4.032 


7.253 


7.696 






Ge 


-0.152 


-0.076 


-0.071 


-0.255 


-0.466 


-0.485 






Gm 


6.213 


5.650 


5.701 


6.509 


5.848 


5.469 






EMR 


-2.45% 


-1.35% 


-1.24% 


-3.92% 


-7.97% 


-8.87% 


N(1440) 


Pll 


h 


0.213 


0.304 


0.303 


0.272 


0.302 


0.300 






9" 


0.255 


-0.269 


-0.247 


0.255 


-0.164 


0.017 








-0.125 


0.273 


0.234 


-0.125 


0.096 


-0.128 


N(1520) 


Di3 


h 


0.560 


0.567 


0.366 


0.560 


0.567 


0.360 






9^1 


-5.753 


-4.848 


-5.607 


-5.498 


-0.580 


-2.348 






9l 


1.217 


2.829 


1.982 


0.301 


-1.503 


0.105 






V 

92 


-0.861 


-0.645 


-0.520 


-0.920 


-0.986 


-0.691 






9^ 


1.462 


0.960 


0.979 


1.674 


2.731 


2.174 


N(1535) 


Sii 


h 


0.132 


0.079 


0.078 


0.083 


0.078 


0.079 






9^ 


0.219 


0.078 


0.028 


0.435 


0.230 


0.084 








-0.102 


-0.127 


-0.080 


-0.164 


-0.195 


-0.129 


A(1620) 


S31 


h 


0.133 


0.159 


0.155 


0.126 


0.159 


0.155 






9 


-0.154 


-0.324 


-0.308 


-0.063 


0.008 


0.044 


N(1650) 


Sii 


h 


0.102 


0.132 


0.107 


0.110 


0.134 


0.107 






9" 


0.113 


-0.167 


0.025 


0.117 


0.074 


0.127 








0.018 


0.411 


0.324 


0.019 


0.281 


0.056 


A(1700) 


D33 


h 


0.285 


0.149 


0.528 


0.285 


0.149 


0.528 






91 


-3.513 


0.663 


-11.875 


-3.996 


-19.642 


-26.531 






92 


1.871 


0.548 


-2.392 


2.000 


3.701 


7.293 



reader with all the parameters of the model. Table 2 shows masses and widths, 
the cutoff A, as well as the vector meson parameters Kp, -F^atat, and K^^, for 
each set. Table 3 provides all the coupling constants of the resonances as well 
as the the E2/M1 Ratio (EMR) of the A(1232) resonance. Table 4 contains 
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Table 4 



Helicity amplitudes in GeV for the different sets. 









#1 


#2 


#3 


#4 


#5 


#6 


A(1232) 


P33 


^1/2 


-0.129 


-0.123 


-0.123 


-0.129 


-0.101 


-0.090 






A A 

"^3/2 


-0.247 


-0.225 


-0.224 


-0.263 


-0.248 


-0.231 


N(1440) 


Pii 


^1/2 


-0.061 


0.064 


0.058 


-0.061 


0.039 


-0.004 






An 

^1/2 


0.030 


-0.065 


-0.055 


0.030 


-0.023 


0.030 


N(1520) 


Di3 


^1/2 


-0.020 


-0.020 


-0.034 


-0.015 


0.027 


0.006 






An 
"^1/2 


-0.050 


-0.013 


-0.022 


-0.068 


-0.121 


-0.092 






"^3/2 


0.161 


0.129 


0.136 


0.161 


0.095 


0.092 






An 

^3/2 


-0.128 


-0.118 


-0.107 


-0.128 


-0.190 


-0.163 


N(1535) 


Sii 


^1/2 


0.060 


0.022 


0.008 


0.119 


0.065 


0.024 






4" 


-0.02(S 


-0.036 


-0.023 


-0.015 


-0.055 


-0.037 


A(1620) 


S31 


"^1/2 


0.038 


0.081 


0.077 


0.016 


-0.002 


-0.011 


N(1650) 


Sii 


"^1/2 


0.037 


-0.054 


0.008 


0.037 


0.024 


0.041 






"^1/2 


0.006 


0.133 


0.105 


0.006 


0.091 


0.018 


A(1700) 


D33 


"^1/2 


0.109 


0.015 


0.222 


0.119 


0.406 


0.573 






^3/2 


0.063 


0.055 


0.057 


0.063 


-0.156 


0.006 



the helicity amphtudes of the resonances which can be compared to those in 
other references such as [11,12,13,30]. 



4-2 Multipole analysis 

As has been previously explained, in order to determine the parameters of the 
resonances and the cutoff we have used the data for electromagnetic multi- 
poles. In this section we discuss the results obtained for multipoles as well as 
the quahty of the different fits. 

Lagrangian models like the one presented here arc more complicated than 
Breit-Wigner models such as MAID [9]. The latter are simple and describe 
accurately experimental observables but do not provide much information 
about properties of the resonances such as the strength of the couplings. 
Breit-Wigner treatment of resonances can be considered naive because each 
resonance contributes only to the multipole with its same angular momentum 
quantum number In this way there is no background from resonances, which 
is very different from Lagrangian models where, for a given resonance, the 
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direct term contributes only to a single spin-isospin channel, while the crossed 
term contributes to different spin-isospin channels as background, and then 
one resonance does indeed affect to the determination of the parameters of 
other resonances. Contributions from crossed terms to the background cannot 
be neglected and there are resonant contributions to several multipoles. For 
instance, N(1520) contributes to i?2^ and Mfl", as expected for a D13 isobar, 
but also contributes strongly to Aff_^. Thus, the background of Breit-Wigner 
models is much simpler because it only has contributions from Born terms 
and vector mesons (p and cu). 

Figures 4, 5, and 6 show the comparison of the six different sets of Table 2 to 
experimental data from SAID database [4]. Without FSI, at low energies, we 
get nice fits to some of the multipoles: , -^2-^, and -Ef-- With increasing 
energy there is a breakdown of the model which calls for further improve- 
ments. The major ingredient that lacks the model is FSI, which we introduce 
phenomenologically as described in section 4.1.3. Indeed, the fits are greatly 
improved - specially the fits of the imaginary parts of the multipoles - when 
FSI are inchided, as it stems from the comparison of the (Table 2). The 
experimental data are quite well reproduced by theory with better quality for 
the low energy region than for the high energy (900 MeV and further), where 
some of the fits start to diverge (i.e. ImMf^ and IuiEq^). In this section we 
focus on fits that include FSI, except in the case in which comparison with 
non-FSI sets provides relevant information. 

Despite of the difference between SP and Vrana et al. masses and widths, the 
curves that we obtain for sets #2 and #3 are very close to each other (so 
are their x^, see Table 2), sometimes undistinguishable, except for some high 
order multipoles as ImMf.. Curves from set #1 do not reproduce data as well 
as #2 and #3 do and the is almost twice as large due to the additional 
restrictions in the values of the parameters. 

If we go through the multipoles in detail, it is convenient to start with Ei^^ 
and (both in Fig. 4) which provide information about the most impor- 

tant low-lying nucleon resonance, the A(1232). These multipoles are of great 
interest at present and a lot of experimental effort has been put in the study of 
the A(1232) in the last years [6,65]. The Mi/^ presents a quite simple struc- 
ture which is very well reproduced by all our sets and it is not affected by 
FSI. That is why all sets are quite similar. Sets ^1, ^A, and ^5 overestimate 
the multipole peaks which will cause an overestimation of the cross section as 
will be seen in section 4.5. The situation is much more complicated for the 
where the FSI are critical, as can be inferred when we compare data 
to sets with and without FSI and check the strong differences among them. 
For these multipoles, data cannot be well reproduced without the inclusion 
of Spsi- When the latter is included the multipoles show a discontinuity at 
= 1.249 GeV due to an abrupt change in SAID phases at that energy. 
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0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

E^(GeV) E^{GeV) 

Fig. 4. Electromagnetic multipoles for the isospin-3/2 channel. Data have been taken 
from Ref. [4]. Photon energy is given in the laboratoy frame. Curves conventions: 
thick dotted set #1; thick solid set #2; thick dashed set #3; thin dotted set #4; 
thin solid set #5; thin dashed set #6. 
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Fig. 6. Electromagnetic multipoles for the isospin-1/2 neutron channel. Same con- 
ventions as in Fig. 4 apply. Data have been taken from Ref. [4] 
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Fig. 7. Examples of various contributions to the multipoles. Left panel shows the 
M^^^ multipole, right panel the Eq^ multipole. Data have been taken from Ref. [4]. 
All the curves have been obtained using set ^2 parameters. Thick dashed: Born 
terms contribution; thin dashed: vector meson contributions; thick dotted: direct 
terms contribution from resonances; thin dotted: crossed terms contribution from 
resonances; thin solid: full calculation without FSI; thick solid: full calculation with 
FSI. 

An important quantity related to A(1232) is the E2/M1 Ratio (EMR) which 
is related to the deformation of the nucleon [6,19,24]. This quantity can be 
defined as the Ge/Gm ratio 



EMR = 



Ge 
Gm 



-- ^^--^)^^ X 100% . (59) 

2(Ma + M)^i + (Ma-M)^2 ^ ' 



We obtain a negative value for this ratio, which according to Ref. [14] corre- 
sponds to an oblate deformation. The values from the most reliable fits (sets 
#2 and #3) are very similar, around —1.3%. This result compares well to some 
other analysis: -1.45% (K-matrix) [11]; -1.42% (ELA) [12]; -1.43% (ELA) 
[14]; -1.3% (dynamical model) [17]; -2.09% (dynamical model) [18]. How- 
ever, it is quite different from the result recently obtained by Pascalutsa and 
Tjon (EMR = 3.8%) [19] within a dynamical model. This ratio is discussed in 
more detail in Ref. [66]. 

The multipoles Mf_ and M{*_ are closely related to the N(1440) resonance. 
If we focus on sets ^2 and ^^3, when Spsi is included the fits look quite 
well except for the real part of the M"_ (second figure of the left panel in 
Fig. 6) where a serious discrepancy between theory and data is found in the 
0.2 — 0.5 GeV energy range. Also, a rather odd behavior in the Mf_ is found 
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between 0.3 and 0.4 GeV (sec Fig. 5), where no experimental data are available. 
For these multipoles related to N(1440) resonance, background and resonant 
contributions are not well established. As a consequence, the parameters of 
the Pii resonance cannot be well determined. These multipoles also show the 
importance of FSI in the model in order to determine the properties of the 
resonances because of the large discrepancies among fits with and without 
Spsi- However, if we focus on sets #1 and #4, FSI do not seem so important 
if the PDG values are used. Actually, set #4 provides better results than set 
#1 except for the high energy region of ReMf_ and ReM{*_. More research on 
the properties of this resonance (and of its role in nuclear medium) has to be 
done in forthcoming years [1] . 

Resonance N(1520) contributes mainly to £^2- ^^'^ ^2- due to its angular 
momentum and isospin. It also contributes sizeably to other multipoles. The 
s-channel contributes to Mf^_ and its crossed term to ImM^^^ as background. 
It also has small contributions to the background of other multipoles. Consid- 
ering set 7^2 and multipoles -Efl" and Affl", the agreement is excellent except 
where there are few experimental data. Set #3 overestimates the peak of the 
resonance in the multipoles and so will do for the cross section. 

^3/2,p,n j-Q^|^jpQ]gg gg|^ contributions from Born terms and vector mesons mainly. 
Resonances N(1535), A(1620), and N(1650) only contribute in the high energy 
region, but in that region they acquire great importance defining the shape of 
the multipoles. For example, the cusp peak that shows up in Im£^Q^ (Figs. 5 
and 6) is due to the structure of the phenomenological width - Eq. (44) - and 
to the inclusion of the partial decay width F^/F in N(1535) resonance. Multi- 
poles i?o^ are well reproduced by sets 7^2 and ^^3, except in the high energy 

region of Imi^Q^. The multipole Eq\_ (Fig. 4) is not so well reproduced in the 
intermediate energy region (0.4 — 0.8 GeV), with an overestimation of the real 
part and an underestimation of the imaginary part. This indicates that the 
prediction of the model is correct for the absolute values of the multipoles and 
that there may be a problem with the phases. 

Only one resonance remains, A(1700), which is associated mainly to multi- 
poles £'2^^ and mI^_^ . As one can see in Fig. 4, when enough data points are 
available the fits are good, yet the large ambiguities in the mass and width 
of this resonance make somewhat unreliable the determination of its coupling 
constants and its contribution to the observables (see Table 3). Further re- 
search on the properties of this resonance is necessary. 

In figure 7 we show two examples of the various contributions to the multipoles 
using the coupling constants of set ^2. It is clear that, without FSI, the 
Born terms and vector mesons do not contribute to the imaginary part of the 
multipoles and represent a background — it has to be noticed that when the 
FSI are included, they do contribute to both real and imaginary parts of the 
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Table 5 



Reduced cross section at threshold in /xb/sr. Experimental data have been 

taken from reference [12]. 



Sets 




#1 


#2 


#3 


#4 


#5 


#6 


Experiment 




pn 


0.0984 


0998 


0998 


0.0949 


0.1023 


0.1020 


094 ± 017 






0.0046 


0.0045 


0.0045 


0.0049 


0.0044 


0.0044 




7n — 




18.92 


18.93 


18.93 


18.91 


18.95 


18.95 


20.4 ±0.7 


















20.0 ±0.3 


















19.7 ± 1.4 






14.51 


14.50 


14.50 


14.52 


14.48 


14.49 


15.4 ±0.5 


















15.6 ±0.5 



multipolc. Left panel shows the multipole M^^ , whose main contribution is 
the A (1232). In this multipole, FSI are not important and curves with and 
without SAID phases differ little. Thus, the phenomenological width included 
is enough to describe accurately the multipole and its structure is quite simple. 
However, the situation is different for the multipole £"0+ which presents a more 
complex structure because its dominant contribution comes from Born terms 
and vector mesons. In the absence of FSI the imaginary part of this multipole 
is practically zero up to 0.8 GeV. While inclusion of FSI makes Born and vector 
mesons contribute to the imaginary part too, improving agreement with data. 

We have not considered spin-5/2 resonances in the model. This will be required 
in order to extend the model to multipoles of higher angular momentum. For 
the energy range considered here, their contribution is expected not to be 
important, although their contribution to the background could improve the 
agreement with data. 

4-3 Results at threshold energy 

Special attention has to be paid to the behavior of the model at low/threshold 
energy, because cross sections and multipoles are well predicted by Low Energy 
Theorems (LET) [67] and Chiral Perturbation Theory (ChPT) [22,23]. Owing 
to the change in the spin-3/2 coupling scheme, the threshold energy results 
change substantially when compared to previous works. In particular, in Ref. 
[12], using the off-shell formalism, it was found that the contributions from 
resonances, direct and crossed terms, were of great importance in order to 
explain the reduced cross section at threshold and the low energy behavior 
of the cross section. These contributions were specially important in neutral 
processes, mainly because of A(1232), and were almost the total contribution 
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to the nn^ production channel [68]. However, in the present calculation we 
obtain a zero contribution to the reduced cross section at threshold from both 
direct and crossed resonance terms. The reason for such a change is the spin- 
3/2 coupling scheme used in the present article, which has no spurious spin- 1/2 
sector. The reduced differential cross section at threshold is proportional to 
the Eo-\- multipole [23] which is a spin-1/2 multipole. Thus, at threshold, any 
contribution of the direct channel from spin-3/2 resonances is a contamination 
which unveils a pathology in the model. This is the case of models based upon 
the traditional spin-3/2 formalism explained in section 3.4.1, as the one used in 
Ref. [12]. This result is independent on the phenomenology of the decay width 
and on the form factors. Therefore, we conclude that, at threshold, only Born 
terms (Fig. 1) and vector mesons (Fig. 2, diagram (e)) contribute, as the spin- 
1/2 resonances are at much higher energy. In Table 5 we present results for 
the reduced cross section at threshold for the four different processes and we 
find a good agreement with experimental values. Tiny differencies are found 
with various parameter sets (#1 to #6) that use different cutoff A and small 
variations in the vector meson parameters. 



4-4 Differential cross sections and asymmetries 



In this section we show results for the differential cross sections together with 
results for five asymmetries: the recoil nucleon polarization, P; the polarized 
target asymmetry, T; the polarized beam asymmetry, S; and the double polar- 
ization parameters G and H. Details on the definition of these quantities can 
be found in Refs. [8,34,69]. The asymmetries are of great interest in the search 
of missing resonances which do not show up so clearly in other observables 
[70]. The formulae which relate the amplitudes with the asymmetries will be 
presented in forthcoming paragraphs. We provide the reader with a wide sam- 
ple of figures in order to have a broad outlook of the model compared to data 
whenever available. 

The FSI treatment described in section 4.1.3 has been applied only to the 
7P 7r°p process. For the other three pion production processes, no FSI 
phases have been included because we have no means to determine them from 
the data available. We calculate the observables for these processes for the six 
sets of coupling constants obtained by fitting 7p — > 7r°p multipoles, given in 
Tables 2 and 3. Thus, these calculations have no adjustable parameters. As 
we shall see in what follows an overall good agreement with data has been 
found. As the energy increases, difi^erences among the curves obtained with 
the different sets of parameters show up more Data favour the sets of coupling 
constants obtained using FSI. 
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Fig. 8. Differential cross section in ^b/sr of the ir^p reaction for different 

photon energies in tlic laboratory frame. Pion scattering angle d in the center of 
mass reference system. Experimental data have been taken from reference [4] and 
are within the range Ej ± 1 MeV. Curve conventions: thick dotted set #1; thick 
solid set #2; thick dashed set #3; thin dotted set #4; thin solid set #5; thin dashed 
set #6. 
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Fig. 9. Same as in Fig. 8. Experimental data are within the range Ej±3 MeV. 
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Fig. 10. Recoil nucleon polarization of the 7p — ^ 7r°p. Photon energy in the labora- 
tory frame. Pion angle in the center of mass reference system. Experimental data 
are within the range E^±3 MeV. Conventions for the curves are as in Fig. 8. 
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Fig. 11. Polarized target asymmetry of the 7p Tr^p reaction. Experimental data 
are within the range zt 3 MeV. Same conventions as in Fig. 10. 
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Fig. 12. Photon beam asymmetry of the 7p n^p reaction. Experimental data are 
within the range Ey ±3 MeV. Same conventions as in Fig. 10. 
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Fig. 13. G asymmetry of the n^p reaction. Same conventions as in Fig. 10. 

Let us first consider the process jp —>■ n^p, for which the experimental database 
has been largely increased in the last ten years mainly thanks to the experi- 
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Fig. 14. H asymmetry of the 7p n^p reaction. Same conventions as in Fig. 10. 

mental programs developed at Mainz (MAMI) and Brookhaven (LEGS). For 
this process the amount of experimental information is much larger than for 
any other pion photoproduction process. Even so, the database on asymme- 
tries is not yet large enough and more measurements are needed in order to fill 
in the existing gaps. Figures 8 and 9 show theoretical curves for the differential 
cross sections compared to experimental data. Differential cross sections have 
been calculated using equations from section 2 and amplitudes from appendix 
A. 

Because of parity, among the eight helicity amplitudes, only four of them are 
independent 



Hi =.4.1/2,-1/2,1 = ~v4_i/2,l/2,-l , 
-f^2 = ^-1/2,-1/2,1 = -4,1/2, 1/2,-1 , 
-f^3 = -4,1/2,1/2,1 = -4_i/2, -1/2,-1 , 
-f^4 = ^-1/2,1/2,1 = ~^l/2,-l/2,-l • 



(60) 
(61) 
(62) 
(63) 



In terms of these four independent helicity amplitudes (see section 2 for 
v4Ai,A2,A-y definition), it is possible to define all the physical observables [2]. 
In particular, the five asymmetries previously mentioned. 

Focusing on sets with Spsi phases, the fits are qualitatively good in the whole 
energy region, and even quantitatively so in the range 250 — 400 MeV. Asym- 
metries are well predicted in almost the whole energy range. 
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Fig. 15. Differential cross section of the 7n — > n^n reaction. Experimental data are 
within the range it 5 MeV. Same conventions as in Fig. 8. 
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Fig. 19. G asymmetry of the 7n — > 'K^n reaction. Same conventions as in Fig. 10. 
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Fig. 20. H asymmetry of the 777. — > vr^n reaction. Same conventions as in Fig. 10. 
In Fig. 10 we provide recoil nucleon polarization asymmetries (P) defined by 

1 k 



ai9)P 



GAn^s* E* 



Im 



(64) 



where the bar over the helicity amplitudes Hj stands for complex conjugate 
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Fig. 21. Differential cross section of the 7p n^n reaction. Experimental data are 
within the range it 5 MeV. Same conventions as in Fig. 8. 
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Fig. 22. Recoil nucleon polarization of the ■yp Tr~^n reaction. Experimental data 
are within the range E'-y it 3 MeV. Same conventions as in Fig. 10. 

and (J {6) for the differential cross section given by Eq. (5). Up to 600 MeV, 
data are well reproduced by sets with FSI. Above this energy, data are repro- 
duced qualitatively but not quantitatively. 
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In Fig. 11 we present polarized target asymmetry (T) given by equation 



a{e)T 



647r2s* E* 



Im 



H2H1 + H4HS 



(65) 



Up to 400 MeV the six curves are very similar. For 500 and 580 MeV the sets 
with phases provide good results and the sets without phases do not. The high 
energy region (700 and 800 MeV) is not well reproduced in general. 

Polarized beam asymmetry (E) is well predicted in the whole energy range by 
sets with FSI (Fig. 12). Even sets without FSI provide good results except in 
the very high energy region (800 MeV). Helicity amplitudes are related to E 
through 

1 k* 



-Re 



H2H3 — H1H4 



(66) 



In short, compared to data, good agreement is obtained for energies below 800 
MeV. Beyond that energy some observables (v.g. E) are also reasonably well 
described 

In the energy region considered here there are no experimental data on the 
other two asymmetries G and H. These asymmetries are expressed in terms 
of helicity amplitudes by means of the following equations 



k* 



647r2s* E* 



Im 



H^H^ + H1H4 



(67) 



a{e)H 



6477^5* E* 



Im 



H2H4 + H3H1 



(68) 



We have also calculated these asymmetries and our results are presented in 
Figs. 13 and 14. 



4.4.2 -fn TT^n 

The situation for the ■jn —>■ 7r°n process is quite different from the previous 
case. The amount of experimental information is very small: No asymmetry 
data are available and the differential cross section data are scant. In figure 

15 we show differential cross sections and in figures 16, 17, 18, 19, and 20 the 
predicted asymmetries (P, T, S, G, and H respectively) obtained with sets of 
parameters. There is a reasonable agreement with data, and sets #1 and #4 
(PDG values) provide the best results globally. 
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Fig. 23. Polarized target asymmetry of the ^ Tr+n reaction. Experimental data 
are within the range ±A MeV. Same conventions as in Fig. 10. 
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Fig. 24. Photon beam asymmetry of the 'jp — > Tr~^n reaction. Experimental data are 
within the range ±A MeV. Same conventions as in Fig. 10. 
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Fig. 25. G asymmetry of the jp Tr~^n reaction. Experimental data are within the 
range E^±3 MeV. Same conventions as in Fig. 10. 



4.4-3 Charged pion production 

In the next paragraphs we go in detail through the predicted differential cross 
sections and asymmetries for charged pion processes, and compare them to 
available data (Figs. 21-32). 

7P Ti^n differential cross sections (Fig. 21) are well predicted by the model 
in the whole energy range by all parameter sets. In the high energy regime (two 
last figures of the panels) differential cross sections are not well predicted by 
any of the parameters sets in the forward scattering region, with the exception 
of set #1 (PDG with Spsi) which provides an impressively good agreement. 
For the P asymmetry (Fig. 22) all curves arc alike and reproduce data correctly 
up to 400 MeV. As the energy is increased, sets #1 and #4 (PDG values) 
provide the best results. The T asymmetry is quahtatively well predicted, but 
quantitative agreement is only achieved up to 500 MeV (Fig. 23). Sets with 
and without FSI provide a good agreement with data for the E asymmetry 
(Fig. 24). Only in the last figure of the panel (700 MeV) wc observe different 
qualitative behaviors from one set of constants to another. Data are scant and 
not reliable for G and H asymmetries. As in previous asymmetries, in the 
low energy regime all the curves are alike, but as energy is increased their 
predictions become quite different. 

Differential cross section data for the reaction 7p — > T{~p are well predicted 
by the sets with FSI (#1, #2, and #3). All the curves are similar for the P 
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Fig. 26. H asymmetry of the 7p ir'^n reaction. Experimental data are within the 
range it 3 MeV. Same conventions as in Fig. 10. 
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Fig. 27. Differential cross section of the 7J) — > tt p reaction. Experimental data are 
within the range E'-y =b 4 MeV. Same conventions as in Fig. 8. 
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Fig. 28. Recoil nucleon polarization of the 7p — > tt p reaction. Experimental data 
are within the range ±1 MeV. Same conventions as in Fig. 10. 

asymmetry (Fig. 28) and are close to data. Overall agreement is good for the T 
asymmetry (Fig. 29). This agreement becomes excellent for the highest energy 
{E^ — 802 MeV) if we consider only curves #2 and #3. E asymmetry (Fig. 
30) is very well predicted by curves #2 and ^^3 in the whole energy range. All 
predictions are qualitatively quite similar for the G and H asymmetries (Figs. 
31 and 32) except for = 800 MeV, where large differences are found. 

The model works quite well for processes with charged pion. So, this is remark- 
able if we take into account that no Spsi have been included, and indicates 
that FSI are not as important in the studied energy region, for charged pi- 
ous as they are for neutral pion channels. Quantitatively, the model provides 
satisfactory results nearly in the whole energy range and in almost every ob- 
servable. Even in the cases where good quantitative result is not achieved, at 
least the qualitative behavior of data is well reproduced (i.e. Figs. 21 and 23). 



4-5 Cross sections 



Finally, in Fig. 33 we present results for the total cross sections compared 
to available experimental data. The two upper figures show the total cross 
section for charged pion channels and the two lower are for neutral pion pho- 
toproduction. 

It is interesting to notice how some of the observed effects in the multipoles 
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Fig. 29. Polarized target asymmetry of the -fp ^ n p reaction. Experimental data 
are within the range ±5 MeV. Same conventions as in Fig. 10. 
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Fig. 30. Photon beam asymmetry of the 'jp ^ n p reaction. Experimental data are 
within the range ±1 MeV. Same conventions as in Fig. 10. 
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Fig. 31. G asymmetry of the 7p — > tt p reaction. Same conventions as in Fig. 10. 
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Fig. 32. if asymmetry of the 'jp ^ n p reaction. Same conventions as in Fig. 10. 

show up in the cross sections. For example, sets ^1, ^4, and ^5 overestimate 
the first resonance region due to the overestimation of M'l!^ peak. On the other 
hand, set #4 presents a cusp peak in multipole Imii^Q_,_, that also shows up in 
the cross section, specially so in the ir'^n channel. The high energy behavior 
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E^{GeV) E^{GeV) 

Fig. 33. Total cross section as a function of photon energy in laboratory frame. Same 
curve conventions as in Fig. 8. 

is well regularized. Nevertheless, it has to be considered that we do not take 
into account resonances D15 and F15 which may change the shape of the cross 
section in the second resonance region. 

The low enery behavior of the charged processes is quite well reproduced by 
all the sets of parameters. Actually, curves obtained with coupling constants 
from sets #1 and #2 agree quite well with data in almost the whole energy 
range. Other sets do not provide good results: sets 7^5 and #6 overestimate 
greatly the second resonance region for Ti~p channel, and set #4 does the same 
in TT+n channel. Overestimation of the second resonance region by set #3 is 
due to the overestimation of multipoles related to resonance N(1520). 

Concerning the 'jn tt^u channel we found several differences among sets 
either in the region of the first or in the region of the second resonance. As no 
data are available for 7r°n total cross sections, we rely on results on differential 
cross section to infere that: Up to 400 MeV, sets #2 and #3 may provide a 
good estimation of the total cross section; and that beyond that energy, there 
may be probably an underestimation of the total cross section. 

In summary, we conclude that set #2 is the most reliable one because it 
provides the best results when all data are considered globally. We weigh that 

this is so, regardless of the fact that other sets may provide better fits to 
individual cases. For instance, set 7^6 provides the best fit to vr^p total cross 
section and set #1 is very good for charged pion channels. As a matter of fact, 
set #2 has the lowest for the electromagnetic multipoles. With this set the 
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only deviations from experimental data in the total cross section are the slight 
underestimation of Tr+n and n^p processes beyond 400 MeV. 



5 SummEiry and final remcirks 

We have elaborated on a pion photoproduction model which is based on an 
Effective Lagrangian Approach (ELA) and is guided by Weinberg's theorem, 
fulfilling chiral symmetry, gauge invariance, and crossing symmetry. We have 
included Born terms, p and u mesons exchange, and seven nucleon resonances: 
A(1232), N(1440), N(1520), N(1535), A(1620), N(1650), and A(1700). Under 
these premises, the model is independent of the underlying subnuclear physics 
(quarks, gluons), which is embedded in the parameters of the model, such as 
coupling constants, masses, and widths. 

With respect to former models along similar lines, this is the first one that 
covers all the well established spin-1/2 and spin-3/2 resonances up to 1.7 GeV 
and, at the same time, fulfills gauge invariance as well as chiral and crossing 
symmetries. Crossing symmetry could not be achieved in previous models 
such as the one of Ref. [12] due, among other things, to pathologies of former 
spin-3/2 Lagrangians. This problem is fixed in the present work by: 

(a) the use of a spin-3/2 Lagrangian due to Pascalutsa that contains no 
spurious spin-1/2 components in the direct channel, 

(b) the use of consistent, energy dependent, strong couplings and widths, as 
well as form factors. 

One of the goals of this paper is to estabhsh a rehable set of parameters for 
the model. In addition to the cutoff A - which is related to short-distance 
effects and can be considered as the only free parameter of the model - we 
adjust electromagnetic coupling constants of the nucleon resonances within 
the usually accepted ranges. The determination of the parameters has been 
performed by fit to the experimental 7p — > 7r°p multipoles, through a min- 
imization procedure. In the minimization we have considered three different 
sets of masses and widths: 

(a) Masses and widths taken from PDG with electromagnetic coupling con- 
stants within the PDG error bars. 

(b) Masses and widths taken from the multichannel analysis of Vrana et al. 
[60] with electromagnetic coupling constants considered as free parame- 
ters. 

(c) Masses and widths obtained by means of a speed plot calculation with 
the electromagnetic coupling constants considered as free parameters. 
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On the other hand wc have considered final state interactions (FSI) phe- 
nomenologicaUy by adding an extra phase to the 7p n^p muhipoles, in 
order to match the current energy dependent solution of SAID phases [4]. In 
all, we have derived six sets of parameters, one with and one without FSI for 
each of the above mentioned sets of masses and widths: 

(a) Sets #1 and #4. 

(b) Sets #2 and #5. 

(c) Sets #3 and #6. 

Electromagnetic multipoles for 7p vr^p are globally well reproduced by 
sets #1, #2, and #3 that include FSI. The fits without FSI (#4, #5, and 
7^6) are also good in the low energy regime. Other experimental observables 
are surveyed such as differential cross section, asymmetries, and total cross 
sections. At threshold we find good agreement with experimental data. In 
our model almost all the contribution at threshold comes from Born terms at 
variance with results in Ref. [12]. 

For charged pion production, where we have no adjustable parameters, the 
agreement is remarkably good for almost all the observables. We note that 
FSI phases obtained for the 'jp Tr^p process are not applicable to charged 
pion production. Thus, no FSI phases have been included in these charged 
pion photoproduction calculations. The fact that we get good agreement with 
data indicates that FSI are small in 7p — > ir'^n and 771 — > 7r~p. 

Although all the parameter sets are reasonable, we favor set #2 because of its 
lowest to the multipolc data and its better agreement with the total cross 
sections for all processes. Set #3, which also has a low x^, is very similar to 
set #2 and also yields similar helicity amplitudes for all the resonances except 
for A(1700). This resonance is poorly known and more precise information 
would be necessary. A better experimental knowledge of multipole would 
improve the determination of the properties of A(1700) resonance. Similarly 
better knowledge of the Mf_ multipole would help to establish more reliably 
properties of N(1440). 

For the future, it would be interesting to analyze contributions from spin-5/2 
nucleon resonances. Although they are not essential to the multipoles con- 
sidered here, they may contribute to the background and their effect can be 
sizeable in the second resonance region of the total cross section and asym- 
metries. The incorporation of the spin-5/2 resonances will require to take 
into account higher order multipoles in the analysis. The inclusion of other 
resonances not considered here (three stars in PDG and misssing resonances) 
could also improve the fits in some energy regions, but it is difficult to perform 
a reliable determination of the parameters without the aid of other physical 
processes where their contribution would be more sizeable. 



49 



The results obtained here are encouraging and estimulatc the apphcation of 
this model to other processes such as pion electroproduction, two pion pro- 
duction in nucleons and nuclei, as well as electro- and photo- production of 
other mesons. 
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A Invariant amplitudes 



In this appendix we show all the invariant amplitudes needed for the calcula- 
tions in the isospin decomposition and the notation for kinematics of section 2. 
We name v to the four momentum of the exchanged particle in each diagram. 



A.l Born amplitudes 

• s-channel (diagram (a) in Fig. 1) 



lO. 



u{p) 



u-channel (diagram (6) in Fig. 1) 



A?, 



2m7r 



u{p') 



S ^2 



2M 



a: = -A- = K (i^f - Fi F^) 
t-channel (diagram (c) in Fig. 1) 

-u{p') ^\ ' ^^^ l^lbuip) 



^- ^ _^pvUN ^^j^,^A- (v + k) 



t — mi 



(A.l) 
(A.2) 



(A.3) 
(A.4) 



(A.5) 
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KroU-Rudermann (contact) term (diagram (d) in Fig. 1) 



(A.6) 



A. 2 Vector meson amplitudes 
• p meson (diagram (/) in Fig. 2) 



t — rrip 



uj meson (diagram (/) in Fig. 2) 



At = -ie """-^^ """ ^(pO 



t - ml 



7. + ^7a/3^^^ 



7a + ^7a^^^ 



k(p) (A.7) 



u{p) {A.i 



A. 3 Sii resonance amplitudes 
• s-channel 

egsh 



4° 



^^^u{p')^G{v)A^^,,q'^^,u{p) 



(A.9) 
(A.IO) 



u-channel 



4° 



^^^^ ^(p')^'^7M^g^75G(^;)|^«(p) 



(A.ll) 

(A.12) 



^..^ 5'3i resonance amplitudes 
• s-channel 

^tsa, = -24:83, = ^^«(pOl^G'(^)A^7M.g'^75^^(p) 



u-channel 



A+ — OA- ^ 

^w,S3i ^^u.Sai 3 Mf 



(A.13) 



(A. 14) 
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A. 5 Pii resonance amplitudes 
• s- channel 



u-channel 



^tpn = -^«,Pn = KpJ9S ^ 9v) 



A. 6 P33 resonance amplitudes 
• s-channel 



^^P33 ~ p,, — 



—ihe 



-S,P33 



UM*M {M* + M) 



X 



u{p) 



u-channel 



j\+ — 2A~ = 

«,P33 «,P33 f^M*M{M* +M) 

X u {p') [igie^yapv^q^A^ + ^27^ {v ■ qA^ - v Aq^) 



A. 7 D33 resonance amplitudes 
• s-channel 

ihe 



A+ — —2A~ — 

«,D33 .,D33 f^MM* (M + M* 



X 



igieujapct>v'^q''A'''-f^ + g2{v- qA^ - v ■ Aq^) u (p) 



u-channel 
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4+ 



2A" 



ieh 



""'^^^ ~ UMM* (M + M*) 
X u {p') [igie^.^apq'^v^A^l^ + g2{v ■ qA^-v Aq^) 



(A.22) 



resonance amplitudes 
• s-channel 

3i/ie 



^0 ^ 

^''^^^ Af^MM* (M + M*) 

i9fe^ap<t>v'^q''A^-i^ + g^{vqAa-v Aq^)] u (p) 



X 



4+ — A~ 



^s,Di3 {91,2 9l,2 



u-channel 



Q 3ieh 

" Af^MM* {M + M*) 

X « {p') [igie^.^am'^v^^'^l^ + 92 {v-qA^-v Aq^) 

^^,013 = ~^w,Di3 = ^°,Di3 (5'f,2 ^ 91,2) 



(A.23) 
(A.24) 



(A.25) 
(A.26) 



B Decay widths 



The coupling constants at the strong vertices are related to the decay widths 
of the resonances. Given the following kinematical definitions 



2M* 



(M*^-M^-miy -4mlM^ 



En = Vk*^ + M2 , 



(B.l) 

(B.2) 
(B.3) 



the decay widths related to the resonance Lagrangians of section 3 are 
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rSsi ^ [E^{Ej, + M) + k 

27rM*/2 2 (£;Ar + M) 



<2l2 



(B.5) 



pD33^_l! ^ I'D ON 

37r/2 M* Etv + M ' ^ ^ ^ 

nf^M* En + M ' ^ ' ^ 



C Electromagnetic Multipoles 



The starting point for multipolc analysis is to define the helicity spinors in the 
cm. reference system for incoming, 




1 


C 



1 





(C.l) 



and outgoing nucleons, 



— sin|, cos|, C'sin|, — C'cos| 
cos|, sin|, C'cos|, C'sin| 



(C.2) 
(C.3) 



where p = E* + M, p' = E'* + M, ( = ^, and C' = y- 

All formulae in this and forthcoming appendices are in cm. reference system 
and with this spinors definition. 

In order to build up the multipoles it is convenient to change the isospin basis 
from A+, A-) to (A^/^, ^A^/^ nA^^^). Both bases are related by means of 
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A^I'^^A^ -A- , 

3 3' 



3 



(C.4) 
(C.5) 

(C.6) 



Defining A = A^ — Ai, initial helicity state along the photon, and = — A2, final 
helicity state along the pion, the spin and isospin projection of amplitudes can 
be written as 

where W = \/s*, j is the spin of the resonance, and o?-^^ {9) are Wigner d- 
functions with the conventions of [30] . The lowest order multipole amplitudes 
are [2,8,71]: 



Ml 





V2 

4 


-"1/2,1/2 + 




\rrI,l/2 
~ [^1/2,1/2 




V2 
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1/2,1/2 




%/2 
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A"" 1/2,1/2 


eL = 


V2 
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1/2,1/2 



/,l/2 
1/2,-1/2. 

o-/,l/2 
-"1/2,-1/2 



W,3/2 
'3/2,-1/2 



8 



/ rr/,3/2 
V^i/2,1 



/2 



if: 



7,3/2 
1/2,-1/2 



1 / rrA3/2 
^ 1^3/2,1/2 ^3, 



3/2 
3/2,-1/2 



(C.8) 
(C.9) 
(C.IO) 

(C.ll) 
(C.12) 
,(C.13) 



D Experimental helicity amplitudes 



In this appendix we present the connection between our amplitudes and the 
helicity amplitudes of the resonances as they are found in Ref. [30] in order to 
relate the coupling constants to the usual partial wave analyses. To perform 
this connection the isospin decomposition (^A^, A^, A^^ is needed instead of 
the one in section 2. Both are related in the following way 
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A^ = f-{A^-A-) , (D.l) 

— ^ (^^ + 2>1~ + 3A°) , (D.2) 
v3 ^ ' 

^" = -L (^A+ + 2A- - 3A°) . (D.3) 



And the helicity amplitudes are given by [12,72] 



7 / 97r k* M* 

where A, yU, j, and {9) have the same meaning as in appendix C; F is the 
total decay width and the pion-nucleon decay width of the resonance as 
defined in appendix B. k* and q* are the pion and the photon momenta in the 
cm. system. We define the kinematical coefficients 



(D.6) 



1 M* a* 

T^- — - fD7) 

4M(M + M*)v^' ' 

to obtain finally the following results: 

• Resonance Sn 

• Resonance S31 

^f/2(S3.) = -^|\/J(M + M-) (D.9) 

• Resonance Pn 

AP%(Pn) = — ^^J^(M + M*) (D.IO) 

• Resonance P33 
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e^y 2 ~ ^^^^ 
^tJ^ [91 + ^92] 



(D.ll) 



(D.12) 



• Resonance D; 
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eT^ [92-^91] 
eT J ^[92 + ^91] 



(D.13) 



(D.14) 



• Resonance D 



13 




(D.15) 



(D.16) 
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